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^ ■ Abstract 

a: 

Q_i| We perform an exploratory study of higher twist contributions to deep inelastic scat- 



tering. We estimate the size of two major sources of higher twist, namely absorptive 
corrections and the vector meson dominance (VMD) contribution. We find that they give 
a sizeable higher twist component of F%. For example at x = 0.01 it is about 8% (17%) 
at Q 2 = 10 GeV 2 (4 GeV 2 ), reaching up to 27% at x = 10" 4 and Q 2 = 4 GeV 2 . At the 
smaller x value the largest contribution comes from absorptive corrections, while at the 
larger values of x the VMD term dominates. 



At large Q 2 the cross section for deep inelastic scattering (DIS) is to a good approximation 
described by just the twist- 2 component of the structure function. That is[] 

° T {l*p) = ~^F T (x,Q 2 ) (1) 

(2) 

with Ft — Ft , where a is the electromagnetic coupling. Indeed in extracting parton distribu- 
tions from DIS data it is commonly assumed that (|l|) is exact and that Ft is given entirely by 
twist-2. To leading order in QCD we have 

F T 2) =J2e 2 g x[q(x,Q 2 ) + q(x,Q 2 )} (2) 
<? 

where q and q are the quark and antiquark distributions and ee q is the charge of the quark. 
For sufficiently small values of Q 2 the higher twist (4,6,...) components of Ft, defined by 

F (4) p (6) 

* = FP + * r + !§ r + ..., (3) 

would be expected to give noticeable contributions to ot- Surprisingly, even though the DIS 
data have become much more precise, the recent global analyses still show no necessity for 
higher twist contributions — despite including data at remarkably low values of Q 2 . Moreover 
parametric fits of F 2 data |l| have found very small values of F 2 at low x. Our objective is to 
explore the role, and to estimate the size, of the higher twist terms. 

^From the point of view of the Wilson Operator Product Expansion (OPE) the higher twist 
terms correspond to operators describing a larger number of partons. Say, for twist-4 we must 
consider operators with four quarks (N\qqqq\N) or four gluons (N\gggg\N), and so on. Strictly 
speaking for each new twist and new operator we have to specify a new input function which 
should be determined by a global fit to the DIS data. Unfortunately the data are not yet precise 

(2) 

enough to determine more than the leading twist, F T , decomposition in terms of partons. 

Rather we will discuss two effects which are expected to give the dominant twist-4 contri- 
butions in the HERA domain. The contributions arise from (a) absorption corrections and (b) 
the vector meson dominance (VMD) contribution. We estimate their size below, but first we 
show their 1/Q 4 behaviour. 

At low x the behaviour of <Jt(j*p) is controlled by the evolution of gluons in the t channel. 
The evolution equation may be written 

( n&\ / n&\ , f 1 dz f Q2 d Q' 2 ®s(Q' 2 ) x x _ /2 . _ , , 

•i dx' fQ' dQ' 2 



hi Q' J 

where the first two terms correspond to the conventional twist-2 DGLAP evolution with gluon- 
gluon splitting function P gg (z), while the last term with the negative sign takes into account 
the higher-twist absorptive corrections^. The value of A is discussed below. The last term 

1 Since we are interested in the possible higher twist effects on parton analyses we concentrate on the cross 
section for the absorption of transversely polarised photons, ctt- 

2 Such a form was first introduced in Ref. ||], and studied further in Ref. |J. 



1 



corresponds to gluon recombination and is shown schematically in Fig. 1, where the four gluon 
twist-4 structure is evident. The extra Q' 2 in the denominator of the absorptive term reflects the 
small probability to find an additional gluon in a small domain of transverse size 1/Q r . Strictly 
speaking the last term contains both twist-2 and twist-4 components coming from evaluating 
the integral at the lower limit Qq and the upper limit Q 2 respectively. The twist-2 component 
may be regarded as an extra input for the twist-2 DGLAP evolution. On the other hand the 
contribution from the upper limit results from the evolution of the four gluon state from Qq to 
Q 2 , and is therefore manifestly twist-4. 

To identify the twist-2 and twist-4 parts of the absorptive correction we simply rewrite the 
last term in (§) in the form 

-A / ••• = -A / ••• + A / ••• . (5) 

JQI JQl JQ 2 

The negative first term combines with the input distribution, xg(x, Qq), in ([|) to give a new 
initial condition for the twist-2 evolution, while the positive second term is the twist-4 com- 
ponent. Thus although the whole absorptive correction is negative, the twist-4 component is 
itself positive. 

A second well-known higher twist contribution to <Jt{i*p) comes from the Vector Meson 
Dominance (VMD) term 

where the sum is over vector mesons V of mass My and where 7y specifies the 7 — V coupling. 
oy(s) is the total Vp cross section at centre-of-mass energy -^i. The <jt{ VMD) component 
dominates crj-^p) at very low Q 2 . The space-time picture of the VMD contribution is com- 
pared with that of the conventional (twist-2) DGLAP contribution in Fig. 2. For the DGLAP 
contribution of diagram (a) the upper quark propagator 1/(Q + k) 2 ~ 1/Q 2 , so that the dis- 
tance between points 1 and 2, Ar 12 ~ l/Q, is very small. Moreover in terms of old-fashioned 
perturbation theory, the invariant mass of the produced qq system M ~ Q. In contrast in 
the small region of phase space corresponding to M 2 ~ M v -C Q 2 the left-hand diagram of 
2(a) may be drawn as in the left-hand diagram of 2(b), where the space-time development 
of the VMD contribution is evident. The photon first creates a light qq pair long before the 
interaction with the target. The presence of the four quark twist-4 structure is evident in the 
schematic right-hand diagram of 2(b), and is reflected by the l/Q 4 behaviour of (||) at large 
Q 2 . In this case, which corresponds to large distances and rather small transverse momentum 
k,T of the quarks, it is better to deal with constituent quarks with masses m u ^ = 350 MeV and 
m s = 500 MeV. 

Absorptive or gluon rescattering effects 

In general the four t channel gluons shown in the lower part of Fig. 1 interact with each 
other. It is convenient to re-organize the perturbative expansion and to consider first the 
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interaction between the gluons in each pair separately and then to consider the interactions 
between the two gluon ladders. The sign and the value of each contribution depends on the 
colour structure. Each pair of gluons may form a singlet ([I]), a symmetric or antisymmetric 
octet (8 s ,8 a ), decuplets (10, 10) or a 27 colour state. The various colour configurations have 
different energy dependences 

a „ ( s ')M0)-2 (?) 

depending on the intercept a(0), where \J~s' is the energy shown in Fig. 1. Information on the 
intercepts comes from the BFKL equation. Now the BFKL kernel K contains two parts — the 
virtual (one-loop) correction, which results in the reggeization of the t channel gluons, and a 
term describing real s channel gluon emission. In general the intercepts may be written in the 
form 

on(0) = 1 + a s [ci(K real ) - (IQrtuai)] (8) 

where the kernels are averaged over the corresponding BFKL eigenf unctions. The virtual 
correction is negative and does not depend on the total colour charge of the gluon pair, while 
the colour factor q for real emission is equal to 3, (|, |) , and -1 for the singlet, octets (8 a , 8 S ), 
decuplets and 27-plet configurations respectively. Due to the Regge bootstrap^ property of the 
BFKL equation, the octet intercepts are as(0) = 1 — the virtual correction cancels the real 
emission part of the kernel. On the other hand from (||) we see that the intercepts of the 10 
and 27 configurations are less than 1 and their contributions therefore decrease with energy. 
Only the singlet- singlet configuration gives an amplitude which grows faster than the twist-2 
contribution as x — > 0, since the singlet intercept «i(0) > 1. 

To estimate the size of the twist-4 term in the gluon evolution equation @j we therefore 
need to evaluate the factor A in the last term for the singlet-singlet configuration. First, it 
contains a colour factor of 9/16 corresponding to the coupling of the four gluon state to the two 
t channel gluons. The best way to obtain this factor is to consider the cross section for heavy 
photon dissociation [j|, |5| and to use the AGK cutting rules ||, which have been justified in 
QCD in refs. 0> §[]■ To explain the remaining content of A it is convenient to write 

A=*-™ (9) 
16 B K J 

The dimensional factor B (which compensates for the extra Q' 2 in the denominator in the last 
term of (^)) comes from the integral over the momentum t = (p — p') 2 transferred through the 
"pomeron" loop (which is indicated by the circular arrow on Fig. 1). In accordance with the 
measurements of heavy photon dissociation at HERA |9j we use 

B = ^P-/ f^dt^ 7.2 GeV" 2 . (10) 
dt I J dt v ' 

Until now we have considered only "elastic" 7*p — > Xp proton interactions. However the 
probability of dissociation of the target proton is not negligible. It is of order 50-70% of the 
3 In other words the octet BFKL amplitude self-consistently reproduces the original gluon trajectory with 

O8(0) = 1. 
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"elastic" interaction [10]. Thus we include a factor C = 1.6 in (||). Finally we include in the 
coefficient A the effect of pomeron- pomeron rescattering (shown in Fig. 3). Such rescattering 
partially fills the rapidity gap and so this contribution is not usually included in the diffractive 
component F 2 D of F2. The contribution increases with energy, but fortunately the dependence 
is rather slow (~ 1 + 0.45(lnl/x) 1//4 ) and so we may include it as a constant K factor. In the 
HERA domain K = 1.63 to 1.76 []11] and so we let K = 1.7 in @. Interestingly our resultant 
value of A is in close agreement with the value (81/16-R 2 ~ 0.2 GeV -2 for R = 1 fm) used in 



phenomenological analyses of absorptive corrections [ 12fl . 

To compute the higher twist AFt contribution to Ft we start with an ordinary (twist-2) 
DGLAP fit to DIS and treat the last term of the evolution equation (f|, |) as a small correction 

^dx' r°°dQ' 2 



xAg(x,Q 2 ) = A f % r^[a s (Q' 2 )x'g(x',Q 

Jx X JQ 2 



/2m2 



(11) 



Now at small x the quark density is dominantly driven by the g — > qq transition. Though a full 
calculation is possible, for our exploratory study it is sufficient to use the approximation due 
to Prytz [O], that is we evaluate 



AF T (x, Q s 



F T {x,Q< 



Ag(2x,Q 2 ) 
g(2x,Q 2 ). 



(12) 



Our pure gluonic estimate of F^ omits the two jet dissociation, 7 — > qq , contribution cor- 
responding to Fig. 4. This graph contributes to the diffractive structure function F 2 D but 
not to the gluon distribution. Rather it contributes to the quark distribution, but with a 
strength which is supressed relative to the corresponding gluon distribution by the colour fac- 
tor (C F /C A ) 2 ~ 1/5. Moreover it gives a higher twist contribution to both Ft and Ft- The 
explicit form of this small correction to F^ has been evaluated in refs. [[14], [151, an d should be 
included in a detailed calculation of the higher twist contribution. 



Vector Meson Dominance contribution 



In the vector meson dominance contribution 
mesons are included 



V=p,LO,(p 



to higher twist, usually only the p, uj and 

M*a v (s) 



7 2 (Q 2 + M 2 Y 



(13) 



We evaluate the pp, up and (ftp cross sections (oy) at centre-of-mass energy ^Js as described in 
the footnote to eq. (p2|). Of course there will be some contribution from the higher mass vector 
meson resonances but there are expectations that they will be suppressed by smaller values of 
l/7y and/or possibly smaller Vp cross sections oy. 

We can also make an alternative estimate of this non- perturbative higher twist contribution 
based on hadron-parton duality. Since for small x the 7* — > qq fluctuations occur over a much 
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longer time scale than the interaction of the qq pair with the target proton, we may use hadron- 
parton duality to write the 7*p cross section in terms of a dispersion relation with respect to 
the invariant qq mass M [[TJJ , 

f°° dM 2 

<1*P) = E I {Q 2 + M2)2 P( s > M 2 )a qq+P (s, M 2 ), (14) 

where cr q q+ p is the cross section for the scattering of the qq system on the proton and where 
the spectral function p represents the density of qq states. We may use perturbative QCD to 
evaluate (|i~4|) . The cross section is given by the probability \A4\ 2 of the 7* — > qq transition 
multiplied by the imaginary part of the forward amplitude describing the gg-proton interaction 

Aqq+p = isa qq+p . (15) 

For transversely polarized photons the amplitude of the 7* — > qq transition is 



KA y^Cj-j - , ^ (e±-fcr)[(l - 2z) A ± 1] y A , + Xm q 5 xy 
Mt = 1 W^W UA(Te±)UA ' = W^r ' (16) 

where the q and q longitudinal momentum fractions and transverse momenta are z, and 
(1 — z), —kx- We use the notation of Ref. [0, which was based on the earlier work of Ref. Jl7| . 
Namely Q 2 and the photon polarization vectors are given by 

Q 2 = z(l-z)Q 2 +m 2 g (17) 
e r = e± = -j=(0,0,l,±i), (18) 

and where A, A' = ±1 according to whether the q, q helicities are ±~. 
In terms of the quark momentum variables we thus obtain 

M7» = E«t| / dz dk 2 T [Z2 + { ' Q 7/^ + ml *cW*. 4) (19) 

where the number of colours A^ c = 3. Before evaluating (0) let us relate this expression to the 
dispersion relation form of fll~4]). We may use 

M 2 = k }±< (20) 
z(l-z) 

to change the integration variable from dk\ to dM 2 . Then (|19]) has the dispersion- like form 

-t(7'p) = ^T.<J * (c JTm' f M2 [*" + (1 " 2)2 1 + 2m «} ' (21) 



<-i 



In comparison to ([14]) we see that (|21|) is a two-dimensional integral. To see the reason for this 
let us consider massless quarks. Then z = |(1 + cos#) where 9 is the angle between the q and 
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the 7* in the qq rest frame. The dz integration is implicit in (14) as the integration over the 
quark angular distribution in the spectral function p. 

We use the additive quark model (AQM) to evaluate a qq+p in fllTf ) which means that each 



quark is assumed to interact with the target proton individually. For forward scattering, (|15l). 
the momentum of the interacting quark is not changed and thus there is no interference — that 
is the initial and final qq states are exactly the same. Now the cross section o" T (7*p) of ( ]T9| ) 
receives contributions from all M 2 up to Q 2 . However to estimate the non- perturbative higher 
twist contribution we must note that the additive quark model result^, 

- 2 

&qq+p — 2°"pp) (22) 

is only valid up to some relatively low mass, M < M . The AQM hypothesis may be justified 
if the separation between the q and the q, Ar ~ 1/A/r, is large in comparison with the qq 
interaction radius R defined by 

c^ el = ttR 2 ~ a p f/9 ~3-5mb. (23) 

As the separation Ar ~ l/k? becomes smaller (that is k? > 1/R ~ 0.5 GeV) the quarks 
start to shadow each other and the cross section o qq decreases. As a result the integration in 
( fi~9| ) is effectively restricted to the region M < Mq, where Mq may be estimated from (|20| ) using 



kx — 0.5 GeV, the constituent quark mass m q ~ 0.35 GeV, and a typical value of z(l — z) ~ 0.2. 
In this way we obtain the upper limit Mq ~ 2 GeV 2 . We impose this limit on the results that 
we present below. 

This part of phase space (that is the domain of small M 2 , small k\, but z ~ 0.5) is not 
that responsible for the leading logarithm behaviour generated by DGLAP evolution. Rather 
the dk\jk\ behaviour of integral comes from the alignment kinematic configuration of small 
z ~ k\jQ 2 (or 1 — z ~ k\jQ 2 \ Then Q 2 of ( |i~7D is of the order of k\. Hence the cross section 
for qq scattering on the proton at large fc^ (that is k\ 3> Aq CD , but k\ Q 2 ) has the QCD 



form a q q +p ~ l/k?, and so (|l^) becomes 

M7» -^J d ~§- (24) 



^From (^) we see that this region of phase space, z ~ k^/Q 2 , corresponds to M 2 ~ Q 2 . 

Thus, in summary, the higher twist contribution comes from the large distance domain, 
Ar ~ 1/kx of small k? and M, but "large" z ~ 0.5, where we have to use a non-perturbative 
(AQM) estimate of cr qq + p - On the other hand DGLAP arises from the perturbative small dis- 
tance regime of large kx and M ~ Q 2 , but small z (or l-z). 



4 To be precise instead of ( |2^ ) we take the <J q q+ P cross sections to be ^\(j(tt + p) + a(ir^p)] for the light qq 
pairs (or the p and ui VMD contributions), and to be a(K + p) + a(K~p) — ^{<t(tt + p) + <r(n~p)] for ss pairs (or 
the 6 contribution). We take the hadronic cross sections from Ref. fLal. 
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Results and Discussion 

Fig. 5 shows the estimates of the higher twist components as a fraction of as a function 
of x for two values of Q 2 . The lower continuous curve corresponds to the effect of gluon 
rescattering and is calculated from fljjD using the MRS(R2) set of partons jl9| . The upper 



two curves correspond to the two alternative ways to calculate the VMD contribution. The 
dotted curve corresponds to AF T calculated from the p, to and contributions to the standard 
VMD formula (0), while the dashed curve is calculated from eq. (0), using the AQM formula 
(22), with the qq invariant mass M restricted to the domain M 2 < Mq = 2 GeV 2 . The two 



estimates of the VMD component are in good agreement with each other. We use constituent 
quark masses (m u d = 350 MeV and m s = 500 MeV) for the VMD estimates. If we were to use 
current quark masses then AFt would be enhanced by about 25%. 

The combined effect of gluon rescattering and the VMD contribution is also shown (by the 
heavy continuous curve) on Fig. 5. To calculate the total effect we use ([34]) and (|12"D for the 
rescattering contribution together with the AQM for the VMD contribution. In the region of 
x = few x 10~ 4 we see that the two higher twist components are of comparable magnitude. 
However the energy dependence of the two contributions is quite different. The 'perturbative' 
rescattering term involves (xg(x)) 2 oc s ' 4 , but actually increases faster than this due to the 
dx' jx' integration. On the other hand the 'non-perturbative' VMD contribution only grows as 
s 08 . Thus as x increases the VMD component becomes the dominant higher twist component. 
The behaviour is evident in Fig. 5. In fact at Q 2 = 4 GeV 2 , for example, by x = 0.1 the higher 
twist contribution (dominated by the VMD component) reaches about 17% of the whole value 
of Ft- On the other hand at very small x the main twist-4 contribution comes from gluon 
rescattering and at Q 2 = 4 GeV 2 and x = 10~ 4 , for example, the total twist-4 exceeds | of the 
Ft value. 

These are very large higher twist contributions indeed, and there is a puzzle why they 
are not evident in fitting to the data. One possibility is that in fitting the data the input 
distribution together with the twist-2 DGLAP evolution is able in a limited kinematical domain, 
say Q 2 ~ 3 — 10 GeV 2 , to mimic the contribution coming from higher twists. If this is the case, 
then we should subtract the higher twist (or at least the VMD contribution which is more or 
less known) from the DIS data before we perform the DGLAP analysis. In other words in the 
fit to the data we should include the VMD contribution (and maybe also gluon screening) as 
well as the conventional DGLAP contributions. An alternative possible resolution of the puzzle 
is that there exists yet another important source of higher twist which enters with a negative 
sign such that it partially cancels out the higher twist component. For instance renormalon 
models have been used to estimate power corrections to deep inelastic scattering [pOfl . 

We conclude that there are theoretical reasons to expect very significant higher twist contri- 
butions to i*2, even for Q 2 as high as 10 GeV 2 . Surprisingly there is almost no evidence for them 
in the dataf]. We have been able to estimate the size of the two major higher twist contributions 

5 In the latest global analysis ElJ partons have been obtained by fitting to data with Q 2 > 2 GeV 2 , and then 
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coming, first, from vector meson dominance and, second, from gluon rescattering in which the 
four-gluon state is composed of two colour singlets. These two higher twist contributions have 
very different energy dependences and there is an increasingly significant resultant higher twist 
effect as x decreases due to the increasing importance of the gluon rescattering term. Of course 
it is possible that F 2 may be subject to other higher twist contributions which partially cancel 
the effects of the two components that we have estimated and which are expected to dominate. 
This exploratory study makes clear that, on the one hand, the higher twist contribution needs 
more detailed theoretical analysis while, on the other hand, it is important to extract as much 
experimental information as possible by fitting the Q 2 dependence of precise data on F 2 (x, Q 2 ) 
in the region of Q 2 < 10 GeV 2 . 
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the analysis was repeated with the omission of data below Q 2 = 10 GeV 2 . The gluon obtained at Q 2 = 10 GeV 2 
from the latter analysis is some 5% higher at small x than the gluon from the former analysis which contained 
low Q 2 data. This could be the first indication that some non-twist-2 contribution is hidden in the data. 
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Figure Captions. 

Fig. 1 An absorptive or gluon rescattering diagram giving rise to a twist-4 contribution. 

Fig. 2 An illustration of the space-time structure of a conventional twist-2 DGLAP contribution 
and a twist-4 vector meson dominance (VMD) contribution. In the latter case the dotted 
lines makes the 4 quark structure evident. 

Fig. 3 An example of an interaction between the colour singlet two-gluon ladders which gives 
rise to the factor K in (|9|). 

Fig. 4 A contribution to twist-4, F T 4 \ associated with the quark, rather than the gluon, distri- 
bution. 

Fig. 5 The higher twist, AFt, contributions compared to Ft at Q 2 = 10 and 4 GeV 2 . Two 
estimates of the VMD fraction are shown; they are essentially coincident at Q 2 = 4 GeV 2 . 
The fraction obtained from the sum of the gluon rescattering and the VMD(AQM) higher 
twist contributions is shown as the heavy continuous curve. 
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